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16.  Abstract 

This  report  is  primarily  concerned  with  the  analysis  of  induced  transient 
'current  and  voltage  pulses  on  buried  shielded  transmission  lines,  due  to 
■earth  conduction  effects  of  rearby  lightning  discharges.  Two  basic  analytical 
[methods  are  presented  in  this  report  to  describe  the  various  kinds  of  coupling 
mechanisms  between  a lightning  discharge  to  ground  and  an  earth-return 
transmission  line.  The  transmission  line  is  assumed  to  be  a long  straight 
horizontal  coaxial  cable  with  an  inner  shield  and  an  outer  armof,  terminated 
on  both  ends  with  typical  communication  equipment  load  impedances.  The 
general  case  is  considered  here,  in  which  the  outermost  conductor  is  not 
necessarily  in  perfect  contact  with  the  conducting  earth  but  has  a contact 
impedance  with  'the  earth,  as  in  cables  with  an  outer  dielectric  covering  for 
corrosion  or  water  protection.  Both  direct  strikes  to  the  cable  via  arcing 
1'ioin  the  terminal  ground  point  of  the  lightning  channel  to  the  outer  conductor 
of  the  coaxial  cable  and  indirect  strikes  to  the  cable  via  conductive  coupling 
mechanisms  through  the  earth  are  considered. 
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FOREWORD 


The  Post-Doctoral  Program  at  Rome  Air  Development  Center  Is  pursued 
via  Project  9567  under  the  direction  of  Dr.  W.  W.  Everett,  Jr.  The  Post- 
Doctoral  Program  is  a cooperative  venture  between  RADC  and  the  participating 
universities:  Syracuse  University  (Department  of  Electrical  and  Computer 

Engineering,  the  U.  S.  Air  Force  Academy  (Department  of  Electrical  Engineering), 
Cornell  University  (School  of  Electrical  Engineering),  Purdue  University 
(School  of  Electrical  Engineering),  University  of  Kentucky  (Department 
of  Electrical  Engineering),  Georgia  Institute  of  Technology  (School  of 
Electrical  Engineering),  Clarkson  College  of  Technology  (Department  of 
Electrical  Engineering),  State  University  of  New  York  at  Buffalo  (Department 
of  Electrical  Engineering),  Florida  Technological  University  (Department 
of  Electrical  Engineering),  Florida  Institute  of  Technology  ( College 
of  Engineering),  Air  Force  Institute  of  Technology  (Department  of  Electrical 
Engineering),  and  the  University  of  Adelaide  (Department  of  Electrical 
Engineering),  in  South  Australia.  The  Post-Doctoral  Program  provides,  via 
contract,  the  opportunity  for  faculty  and  visiting  faculty  at  the  participating 
universities  to  spend  a year  full  time  on  exploratory  development  and 
operational  problem-solving  efforts  with  the  post-doctoral s splitting 
their  time  between  RADC  (or  the  ultimate  customer)  and  the  educational 
i nst i tut  ions . 

This  effort  was  conducted  via  RADC  Job  Order  No.  9567  0006  for  the 
Federal  Aviation  Administration.  Hr.  Fred  Sakate  was  the  FAA  focal  point, 
and  he  participated  closely  in  the  technical  cocitlination  meetings  and 
cable  testing  sessions. 
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CHAPTER  1 
Introduction 


Earth  conduction  problems  are  encountered  in  various  communication 
and  power  system  circuits  in  connection  with  strong  electromagnetic  dis- 
turbances caused  by  nearby  lightning  discharges.  Excessive  interference  in 
conventional  exposed  metallic  communication  lines  is  Indicated  unless 
adequate  protection  measures  are  provided.  This  may  require  extra  electro- 
magnetic shielding  of  certain  important  communication  lines  and  associated 
buildings  housing  sensitive  equipment  and/or  the  installation  of  protective 
devices  on  certain  communication  equipment. 

This  study  is  primarily  concerned  with  the  theoretical  analysis  of 
the  resultant  circuit  disturbances  caused  by  earth  conduction  effects  of 
lightning  discharges.  Two  basic  analytical  models  are  presented  in  this 
study  to  describe  the  various  kinds  of  coupling  mechanisms  between  a 
lightning  discharge  and  an  earth-return  transmission  line. 

The  theoretical  study  described  In  this  report  is  part  of  a larger 
study  program  to  provide  protection  for  communication  electronics  equip- 
ment against  transient  electromagnetic  disturbances.  Current  and  voltage 
pulses  are  induced  in  cables  running  between  buildings  or  equipment 
enclosures.  These  currents  and  voltages  are  then  coupled  into  the  terminal 
equipment.  The  electromagnetic  disturbances  may  be  the  result  of  nearby 
lightning  activities  or  man-made  electromagnetic  pulses.  This  larger 
study,  known  as  the  FAA  Lightning  Protection  Study,  has  been  performed 
by  the  Post-Doctoral  Program  through  several  of  its  member  universities 
for  the  Federal  Aviation  Administration.  The  institutions  include  the 
Air  Force  Institute  of  Technology,  Florida  Institute  of  Technology, 

Georgia  Institute  of  Technology,  and  Purdue  University.  The  Individual 
participants  In  the  FAA  Lightning  Protection  Study  are  listed  In  Appendix  D. 
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1.1  FAA  Lightning  Protection  Study 

Increasing  use  of  solid  state,  integrated  circuit  electronics  in  FAA 
communication  and  control  equipment  means  that  reliance  on  the  over-voltage 
protection  adequate  for  higher  voltage  electron  tube  and  discrete  transistor 
circuitry  would  be  inadequate.  The  over-voltage  protection  of  carbon  blocks 
in  the  several  hundred  volt  range,  and  neon  bulbs,  with  long,  relatively 
high  inductance  leads  in  the  40  - 100  volt  range,  is  not  adequate  for  the 
solid  state  circuits  which  operate  at  lower  voltage  levels  (presently  down 
to  5 vol ts) . 

The  first  phase  of  the  program  is  an  overall  study  and  consists  of 
three  technical  tasks:  (1)  the  determination  of  the  voltage  and  current 

levels  likely  to  be  conducted  to  FAA  equipment;  (2)  the  determination  cf  the 
susceptibility  levels  of  FAA  Instrument  Landing  System  [AN (GRN-27(V) ] ; and 
(3)  the  determination  of  lightning  protective  dev.ces  that  are  available  to 
reduce  the  levels  of  (1)  to  those  permitted  by  (2).  These  three  tasks  have 
been  performed  in  parallel  with  close  interaction  and  are  essentially 
completed ^ ^ ^ ^ . Appendix  0 lists  the  schools  having  primary 
responsibility  for  each  of  the  tasks.  This  report  is  the  result  of  the 
work  done  under  the  first  technical  task.  Only  the  theoretical  foundation 
and  analysis  are  presented  in  this  report;  detailed  numerical  calculations 
are  to  be  reported  in  a companion  report. 

1.2  Lightning  Induced  Transients  on  Burled  Cables 

Numerous  Interference  and  protection  problems  are  encountered  In  the 
development  and  operation  of  extensive  communication  and  power  systems. 

These  problems  are  caused  by  the  Internal  coupling  of  such  systems  with 
each  other  and  by  the  external  presence  of  the  earth  which,  In  some  measure, 
is  involved  as  a return  conductor.  The  earth  also  serves  as  a return  con- 
ductor for  lightning  currents,  which  often  occasion  disturbances  in  commun- 
ication and  power  circuits.  Lightning  disturbances  are  largely  atmospheric 
phenomena  governed  by  the  physical  properties  of  the  air.  However,  the 
behavior  and  effects  of  the  lightning  near  the  surface  of  the  ground  In 
communication  and  power  systems  are  primarily  earth  conduction  problems 
caused  by  the  finite  conductivity  of  the  earth.  Therefore,  problems  arise 
both  In  communication  and  power  system  circuits  concerning  the  protection 
of  transmission  lines  and  associated  equipment  against  interference  and 
possible  breakdown  caused  by  excessive  voltage  or  current  surges  caused  by 
lightning  discharges. 


To  deal  adequately  with  such  problems,  It  is  necessary  to  consider 
theoretical  solutions  to  the  basic  problem  tn  which  the  earth,  as  well  as 
conducting  current  paths,  are  involved  in  the  lightning  discharge.  The 
analysis  of  such  problems  is  inherently  more  complicated  than  the  problem 
of  completely  metallic  circuits  embedded  in  an  insulating  medium,  since  the 
great  extent  of  the  earth  necessitates  the  use  of  electromagnetic  field 
theory,  rather  than  conventional  transmission  line  or  circuit  theory,  In  the 
solution  of  most  aspects  of  the  problem.  It  is  necessary  to  restrict  the 
analysis  to  fairly  simple  fundamental  cases,  in  which  simplified  models  of 
the  earth,  cable,  and  lightning  channel  geometries  are  used,  on  account  of 
the  complexities  that  would  otherwise  arise.  Therefore,  ionization  effects 
caused  by  high  induced  voltages  or  el  ect  i*ol  vt  i c actions  are  not  considered. 

Also,  the  heterogeneous  character  of  the  earth  as  a conductor  and  an  elect- 
rolyte are  not  considered.  Furthermore,  the  extremely  variable  nature  of 
the  lightning  currents  and  voltages  are  not  considered;  however,  typical 
average  values  of  the  lightning  channel  paramo  let's  are  used. 

1 . 3 Overview 

This  report  is  organized  as  follows: 

First,  the  basic  equations  which  govern  the  behavior  of  the  electromagnetic 
disturbances  caused  by  earth  conduction  effects  of  lightning  discharges  are 
listed  In  Chapter  2,  along  with  the  definitions  of  the  scalar,  vector,  and 
Hertz , potent  I al s , which  will  be  used  to  determine  the  form  of  the  solutions. 
Then,  the  fields  due  to  an  electric  dipole  In  free  space  are  given  in  Chapter  3* 
These  free  space  dipole  fields  are  then  generalized  to  determine  the  fields 
of  a vertical  or  horizontally  oriented  dipole  above  a flat  earth.  A knowledge 
of  these  fields  Is  essential  for  determining  the  mutual  coupling  between  the 
dipole-like  lightning  channel  and  a burled  wire.  Next,  the  actual  current 
induced  on  a buried  wire  Is  determined  approximately  (in  Chapter  k)  and 
exactly  (in  Chapter  5)  from  a knowledge  of  the  mutual  coupling  Impedances, 
which  were  determined  previously  in  Chapter  3.  In  both  the  approximate  and 
the  exact  formulations,  the  results  are  presented  in  the  form  of  an  equivalent 
distributed  transmission  line  model  of  the  coupling  phenomena.  In  Chapter  h , 
the  Induced  electric  field  Intensity  along  the  outer  conductor  of  the  cable 
is  also  determined  for  a lightning  stroke  to  ground.  In  Chapter  5,  the 
characteristic  equation  of  the  transmission  line  is  solved,  via  Integral 
transform  techniques,  to  determine  the  characteristic  values  of  the  trans- 
mission line,  e.g.,  the  propagation  constant  and  associated  propagation  modes. 
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The  result Inq  Induced  currant  end  voltage  surges  on  the  outer  armour  of  | 

the  cable  It  thereby  determined  for  a direct  strike  via  arcing  to  the  | 

cable  and  for  an  Indirect  strike  via  conductive  energization  through  the  U 

Hjj 

earth,  Once  the  Induced  current  and  voltage  surges  are  known  on  the  outer  ■’ 

armour  of  the  cable,  the  voltage  and  current  standing  waves  on  the  trans-  | 

mission  line  can  ba  determined  by  finding,  and  Integrating  over  the  length  |] 

of  the  line,  the  Green's  Functions  for  a distributed  voltage  or  current  ij 

source,  as  Is  done  In  Chapter  6.  Next,  the  Induced  current  and  voltage  j; 

surges  on  the  center  conductor  of  the  coaxial  cable,  resulting  from  the  pene- 
tration of  the  lightning  Induced  discharge  through  the  outer  armour  anti  the  ‘ 

‘1 

Inner  shield  of  the  coaxial  cable  are  determined  via  the  Impedance  transfer  : 

functions,  which  are  developed  In  Chapter  7.  Some  of  the  details  of  the 

above  analysts  are  presented  In  the  Appendices. 


CHAPTER  2 
Basic  Equations 


This  chapter  of  the  study  lists  the  basic  equations  which  govern  the 
behavior  of  the  resultant  electromagnetic  disturbances  caused  by  earth 
conduction  effects  of  lightning  discharges.  For  a complete  exposition 
of  the  subject  matter  and  for  other  aspects  of  the  subject  matter  than 
are  of  primary  concern  here,  reference  is  made  to  the  literature  on 
electromagnetic  theory,  transmission  line  and  circuit  theory. 


2.1  Maxwell's  Equations 

We  start  with  the  time  dependent  Maxwell's  equations  for  the  field 
vectors  i,  ti , and 

- vxfcf.t)  - 


vx»(r.t)  - +s|-J(nt> 

V‘3(r,t)  * 0 


(2.1) 


(2.2) 


(2.3) 


V‘tl(r ,t)  = ps  (r  ,t) 


(2.M 


where  p and  J are  respectively  the  source  charge  and  current  densities 
and  J is  the  conduction  current  density.  Obviously,  we  will  be  dealing  with 
four  different  media:  air,  earth,  conducting  wires,  and  insulating  layers. 

Each  medium  is  characterized  by  its  permeability,  permittivity,  and  conductivity: 
ji . , e,,  and  a.,  where  I stands  for  either  a(air),  e(earth),  c (conductor) , or 


I ( insulator) . The  constitutive  relations  tn  each  medtum  are 
t(r  ,t)  * y,#(r,t) 


(2.5) 


$(r,t)  = £j?(r,t) 
)(r,t)  = 0.?(r,t) 


(2.6) 


(2.7) 


It  Is  desirable  to  transform  equations  (2.1)  - (2.7)  to  the  "frequency1' 
domain.  For  this  purpose,  we  define  the  Fourier  transform  pairs  between 
t and  w: 


t(t,  d - 


i 

72tt 


(2.6) 


?(r,u)  » t Cr ,t)e"jwt  dt 

*4»  J-  00 

where  F and  F are  any  of  the  field  variables  or  sources  defined  above 
and  j ■ J-T.  Therefore,  In  the  "frequency"  domain,  the  equations  become 


(2.8*) 


Vx?  = -juS 

(2.1  ') 

Vxfi  = 3 + J + jojD 

(2.2’) 

o 

it 

■J-CO 

t> 

(2. 3') 

v«p  * ps 

(2.4*) 

t-vfi 

(2.5') 

5 - GjE 

(2.6') 

H - 0{t 

(2.7') 

where  the  arguments  r and  w are  dropped,  for  simplicity.  While  the  source 
terms  Ps  and  1 are  arbitrary,  they  are  related  by  the  conservation  of 
charge,  1 .e. 

7-  (J  + J)  = -Jwp  (2.9) 

s s 

Making  use  of  the  conservation  of  charge  (2.9)  and  the  conatltutive  relations 
(2.5* ) - (2.7') » the  flux  densities  0 and  B can  be  eliminated  in  favor  of  the 
field  Intensities  E and  H,  l.e. 


VxE  <*  -JwUj  H 

(2.1”) 

VxH  * +1  we.  E + J 

(2.2") 

1 s 

7.  it  * 0 

(2.3") 

1 

(2.4") 

joie, 


where  the  complex  permeability  and  complex  permittivity  are  defined  by 


~ .ai 

e,  = e,  - j — 


are  complex  quantities. 


2.2  Hertz  Vectors 

Hertz  showed  that  under  ordinary  conditions,  the  field  Intensities  It 
and  H are  derivable  from  a single  vector  function  ir,  known  as  the  Hertz  vector. 
It  can  be  shown,  by  direct  substitution,  that  equations  (2.1")  - (2.V)  are 
satisfied  by  writing 

1 = k2  tt+  V(V-tt)  (2 


H = jojEj  VXTT 

with  the  Hertz  vector  tt  specified  by  the  differential  equation 
V2TT  + k2TT  ® -3  / Owe  j) 


where 


k2  - (3 j + ja?)2  - ioVjCj 


withaiO,  3 > 0. 

It  is  recalled  that  the  general  solution  of  an  inhomogeneous  linear 
differential  equation  consists  of  two  parts:  a homogeneous  solution  and  a 

particular  solution.  One  possible  way  to  express  the  particular  solution  is 


-jkj  |r-r 1 


The  homogeneous  solution  can  be  expressed  In  various  forms.  For  the  problem 

at  hand,  it  is  convenient  to  express  it  in  terms  of  the  elementary  cylindrical 
(4)  -v 

wave  functions  Zn  . More  specifically,  each  component  of  the  Hertz  vector  tt 

can  be  wr I tten  as 


7 


y2  2 

T -kj  and  where  the  dummy  variable  T of  the  Integration 

represents  the  transverse  part  of  the  propagation  constant.  To  Insure 

convergence  at  z-^fc05,  we  choose  Re  Sj  > 0 . 

In  equation  (2.H»)»  the  coefficient  C (t)  is  left  unspecified.  To 

n 

complete  the  description,  it  is  necessary  to  state  the  boundary  conditions: 
the  tangential  components  of  the  field  intensities  1 and  ft  are  continuous 
at  a boundary  free  of  surface  current.  The  boundary  conditions  can  also 
be  expressed  in  terms  of  the  Hertz  vector  tt  and  its  derivatives  and  will  be 
discussed  in  detail  later. 

2.3  Vector  and  Scalar  Potentials 

The  Hertz  vector  is  also  related  to  the  vector  and  scalar  potential 
t,  and  0 via  the  following  relations 

* Pj  (Oj  + jU)£ j ) TT  = jOJEjUjit  (2.15) 

(f)  =»  -V*7r  (2.16) 
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CHAPTER  3 

Fields  Duo  To  An  Electric  Dipole 


3jJ_  An  Electric.  Pipe  la  in  Free  Space 

Consider  an  electric  dipole  ld£  surrounded  by  free  space  (etseQ»  a*0) . 

For  convenience,  the  coordinate  system  is  so  chosen  that  the  dipole  is  located  at 
the  origin  and  oriented  along  the  z-axis.  If  the  medium  is  assumed  to  be  infin- 
itely large,  then  the  only  condition  that  the  Hertz  vector  must  satisfy  is  the 
radiation  condition  at  infinity.  Therefore,  the  Hertz  vector  is  obtained  from 
equation  (2.13)  by  noting  that  Js>=ld£<$(r  )z,  l.e. 


' 2 2 
r + z 


-v  L A . !d£ 

11=1  + 2 JlnSel 


/ 2 2 
/ r + z 


where 


“•Vo- 


it  is  instructive  and  convenient  to  express  equation  (3.1)  as  the  super- 
I ..  . I ; i on  of  elementary  cylindrical  waves.  For  this  purpose,  use  is  made  of 
M-'.f'i.v'.  fe id's  formula ^ ^ * and  we  obtain, 


-►  A.  IdS, 
■n  = + Zl  wT 


...  - /t2 

C;,  = / T - 


J0(rr)e 


' z|5. 


and  i.i'/e  branch  cut  Is  chosen  such  that  - tt/2  < arg  £q  < . A comparison  of 

v : .or  (3.2)  with  equation  (2.14)  shows  that  only  the.  term  independent  of  angle 
<n  T>)  ; used,  as  expected  from  physical  considerations. 

■ Vai-tical  Electric  Delpole  Above  a Flat  Earth 

Mow  consider  two  regions  separated  by  an  Interface  za0.  The  space  above  the 

: 1 1 >• f ...  c:  c*  is  taken  to  be  free  space  (v=y^ , cf=0)  , while  the  medium  below  the 

in<:.>i  late.  is  taken  to  be  earth  (p  -p„,  e , a ).  Our  aim  Is  to  derive  the  expres- 

e 0 e e 

si ens  for  the  fields  generated  by  an  electric  dipole  ld&  oriented  along  the  z-axis 
and  located  at  z=Zq,  as  shown  In  Fig.  la.  ,The  procedure  Is  as  follows.  In  the 
■■*>■■■  of  the  earth,  the  field,  referred  to  as  the  primary  field,  is  due  to  the 
Jit o'e  alone,  and  the  Hertz  potential 


WWWMM* 


Is  given  by  equation  (3.2),  with  the  dipole  shifted  from  z«0  to  z-z0.  In  the 
presence  of  the  earth,  the  primary  field  Is  partially  reflected  by  the  Inter- 
face and  partially  transmitted  Into  the  earth.  It  will  be  demonstrated  later 
that  for  the  geometry  under  consideration,  all  the  boundary  conditions  can 
be  satisfied  by  working  with  only  the  z-component  of  the  Hertz  vector  tr,  and 
only  the  term  Independent  of  angle  (n«0)  In  equation  (2.14)  Is  needed.  Thus, 
for  z > 0,  the  z-component  of  the  Hertz  potential  can  be  written  as 


za 

and  for  z < 0 


. , Idfc 
**  Vnwef, 


,00 


Jo 


te 


+ R (x)e  0 
2* 


JQ(tr)dT 


(3.3) 


ze 


Tzz(t)e 


z? 


e d0(tr)  dx 


(3.^) 


where  R^Ct)  and  “^(t)  are  the  reflection  and  transmission  coefficients  to  be 
2 

determined  and  - jo)pn(o^  + Jwe^) 

\ 


e " ' 0 e " e 
& 


£ =>  /x"  - k 

e e 


tr 


and  the  branch  Is  chosen  such  that— ^ < arg$g  < %• 

To  determine  the  reflection  a.:d  transmission  coefficients  R2z(t)  and 
T (t),  It  Is  necessary  to  appeal  to  the  boundary  conditions.  Since  ffj  “ z it 


zi 


and  it  ^ Is  Independent  of  <p,  we  have,  from  equations  (2.10)  and  (2.11) 


a9  f | ***  3 TT  . n 

? ‘ r-5rtiz  + z <-pr*  k1  *«1> 


H ■ - (dj  + Jwej)  4> 

Thus,  the  boundary  conditions  at  z»0  require  that 


3tt 


zl 


+ Jwee* 


3lT 


ze 

IT 


(3.5) 

(3.6) 

(3.7) 

(3.8) 


>,:V;/y.wUvn ./  f ^ 
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H8C 


for  all  values  of  r.  Since  equations  (3.7)  and  (3.8)  must  be  valid  for  all 
values  of  r,  these  two  equations  can  be  simplified  by  Integrating  with  respect 


to  r , 1 . e . , 


3ir 


za 


3z 


3tt 

ze 

1z“ 


2 2 
k„  it  - k L it 
0 za  e ze 


(3-9) 

(3.10) 


In  equation  (3.10),  the  coefficients  are  expressed  In  a convenient  form  by 
noting  that  p - p r p/s.  By  substituting  equations  (3.3)  and  (3.4)  Into 

3 6 v 

equations  (3.9)  and  (3.10),  we  obtain 


RzZ(x) " h 


k.‘*o + koV 


‘Z0^0 


(3.11) 


Tzz(T)  “ rr 


2ko2T 


“zo^o 


(3.12) 


ke  + k0  *e 


Upon  substituting  equation  (3.11)  Into  equation  (3.3)*  and  making  use  of 
Sommerfeld's  formula,  we  obtain 


za 


+ J dir 


-JV 


•Jk0R 


+ BV 


- A 


(3.13) 


where 


v^r2  + (z  - zQ) 2 


/2  V 


(z  + z0)‘ 


and 


A « 2k, 


r 


° 5o[ke250  + k025e5 


'50(z  + *0)  j (rr)  dr 
e 0 


(3.H) 
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The  field  Intensities  In  the  air  are  determined  by  substituting  equation 
(3.13)  Into  equations  (3.5)  and  (3.6).  More  specifically,  the  cylindrical 
components  of  the  fields,  for  2 > 0,  are 


IdJt  3 7 

+J  TmSq  zt 


(3.15) 


E 

z 


0 


+J 


Idi 


and 


H =0 
r 

U IdU  3 _ 
3?zt 

H «0 
z 


where 


-Jk0R 


-Jk0R' 


- A, 


and 


Z*  " ko  zt 


(3.16) 

(3.17) 


(3.18) 

(3.19) 

(3.20) 

(3.21) 


The  significance  of  the  terms  Z^.  and  are  discussed  In  Appendix  C. 

3.3  Horizontal  Electric  Dipole  above  a Flat  Earth 

Next,  consider  a horizontal  electric  dipole  ld£  oriented  along  the  x~ 
axis,  as  shown  In  Fig.  lb.  If  only  the  x-component  of  the  Hertz  vector 

-V  ■+  * 

Tr,  l.e.,  tr  «=  x tt^,  were  used j the  boundary  conditions  at  the  Interface  z»0 
cannot  be  met  In  fact,  we  have  to  supplement  xir^  with  the  homogeneous 

solution  as  given  in  equation  (2.14)  for  ztt^ . Since  the  boundary  conditions 
are  somewhat  complicated,  It  Is  convenient  to  work  with  Cartesian  coordinates 
first,  and  then  to  transform  the  results  Into  cylindrical  coordinates. 
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a 


S I nee 

TTj  “ ^TTx{  + 2 TT^  (3.22) 

then  the  rectanqular  components  of  the  field  Intensities  1 and  ft  are,  from 
equations  (2.10)  and  (2.11), 


, 3\<i  »S| 

Ex  ■ kt  ’'x, + -£j—  * ssr 

(3.23) 

8Si  32”zi 

C B — ~ + , 

y 9x3y  3y3z 

(3.24) 

„ , 2 . 3\<l  . *\l 

\ ~ kl  uzl  + 3x3z  + 8z2  ' 

(3.25) 

3\l 

Hx-(VJ“l>  3y 

(3.26) 

Sxx,  3xzj 

Hy  “ <VJ“e|)(_5z 5x*' 

(3.27) 

3TX| 

Hz  --(0,+jue,)  8y 

(3-28) 

Thus  the  continuity  of  E , E , H , and  H at  z«0  dictates  that 

x y x y 


3tt  <)ir_  w 

ko  +3T  (-sr + ~5r>  ■ k.  v * ?£  (-sr + ir1 


(3.29) 


a S'"  3*,-  a 

Pf  _ xa  zas  _9_  / xe  ze. 

8yl  3x  3z  ; 3y  1 8x  ^aT; 


(3.30) 


i#- 
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<*waaa«»wigiga^ 


0 Sir 
2 za 


\>  ~§r 


0 a-n 

2 ze 


9y 


(3.30 


« 3tt 
k 2 (—21 
K0  1 3z 


3tt 


„ 3tt 

_ k 2 /.  _xe 

3x  7 e 3z 


za 


) 


3tt_ 


ze 


3x 


) 


(3.32) 


These  equalities  must  be  valid  for  all  values  of  x and  y.  Thus,  by  Integrating 
equations  (3.30)  and  (3.31)  with  respect  to  y,  we  obtain,  at  z»*0, 


3tt  Stt 

xa  . za 

“ST  + “5T 


3tt 


xe 


3tt 


ze 


3x 


3z 


(3.33) 


. 2 .2 
k tt  ■ k it 
0 za  e ze 


(3.3*0 

Making  use  of  these  two  relations,  equations  (3.29)  and  (3.32)  can  be  simplified  to 


,2  .2 
krt  tt  *=  k tt 
0 xa  e xe 


(3.35) 


, 2 3lTxa  _ , 2 9lTxe 

ko  X ■ TT 


(3.36) 


In  the  absence  of  the  flat  earth,  the  primary  field  tt  Is  simply,  from 

Xfl 


equation  (2.13)  and  Sommerfeld's  formula, 


xa 


+j  .JiL 

J A TRUE, 


.co 

J,  "^J°('tr,e 


I z"zg I 


dx 


(3.37) 


In  the  presence  of  the  flat  earth,  there  also  exists  the  zrcomponents  of  the 


Hertz  vector  In  the  air  and  In  the  earth,  l.e.  tt  and  tt  , and  the  x-components 

za  za 


of  the  Hertz  vector  In  the  air  and  In  the  earth,  l.e.  tt  and  tt 

xa  xe 


The  x- 


component  of  the  Hertz  vector  In  the  air  tt  Includes  the  primary  field 

Xfi 


tt  , as  given  In  equation  (3.37).  Before  writing  down  these  terms,  It  Is 

Xa 


necessary  to  examine  their  angular  dependence.  It  Is  noted 


that  it  and  tt  must  satisfy  equations  (3.35)  and  (3.36)  for  z«0  and  for 
Xfi  xe 


all  values  of  0 S r i ",  OS  S tt.  It  Is  also  noted  that  the  primary  part  of  tt  , as 

xa 


given  by  equation  (3.37).  is  Independent  of  <j>.  "Thus,  In  view  of  the  orthogonality 


Jn<fr 


properties  of  eJ  , 0 < <j>  < 2ir,  we  conclude  that  v and  tt  oontaln  only 

xa  xe 
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>;j  jy.g  iiji"h  i ” 


-•*  - 


the  terms  Independent  of  angle  (n«0) , Thus,  for  tt  and  tr  , only  the  n-0 

xa  xe 

term  of  equation  (2.H)  Is  retained,  1a  for  2 > 0 


ld£ 


xa  *4TU0£ 


7T  » 


+Ji; 


"V 


■|z'zol5o  „ , , ^‘•o 


Jq(t r)  dx 


(3.38) 


xe 


1 m r 

J 5moeT  j 

0 J 0 


z£ 

T^Me  * J0(*r)  dr 


(3.39) 


To  determine  the  <f>  dependence  of  ir  and  tt  , equation  (3-33)  Is 

z a z e 

expressed  Ir.  cylindrical  coordinates 

3tt  9tt  Bn 

cos  <j)  -£2-  + - cos 4>  ~5~  + -gf1  (3.33’) 

As  tt  and  n are  Independent  of  <t>  and  equation  (3-33*)  Is  valid  for  all 
xa  xe 

values  of  0 * <p  s 2tr,  tt  and  ir  nust  have  cos <f>- type  dependence.  Thus,  In 

Zd  ZC 

writing  down  the  expressions  for  n and  n , It  Is  only  necessary  to  keep  the 

Z6  Z0 

n«*+ 1 terms,  l.e. 


^za  “ +^Trwe0  ° (Tr)dT^jcos  d> 

» Id2_  zl  -r 

ze  -Hmueg  TX2(T)e  * J,(Tr)dTJcos  <p 

Upon  substitution  of  equations  (3 . 38) - (3 • ^ t ) Into  equations  (3.33')* 

we  obtain: 


(3.40) 


(3M) 

(3.34)-(3.36), 


■r  / \ „ O T 0 0 

-(t> ' 2 r2-  6 

e 


R (t)  - ^ 
xx  C0 


1 - 


2*e  "j 

FZ1 


:*  0*0 
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(3.A2) 


(3.^3) 


»K<« 


rxZ(t)  ’ ‘2  —r — *~z — e 

ke«0  + ko'. 


Thus,  the  Hertz  potential  due  to  a horizontal  dipole  I d 5t  oriented  along  the 
x-axts  and  located  at  z“Zg  can  be  determined  by  combining  equations  (3.42)  - 
(3.45)  with  (3-38)  - (3.41).  It  Is  noted  that  it  can  be  written  In  a rather 
simple  form: 


- ^knR  -Jk0R' 

„ . ..  Id*  e 0 , e 0 A 

^xa  ^ TrfwcT  L R R*  ^ 


']• 


(3.46) 


where  R and  R‘  are  defined  In  Section  3.2  and 


A'  - 2 


f V 


«> 


■Kn(z  + zj 

0 J (xr)  dx  . 
o 


(3.47) 


_ . , I dH  , ' ' 

1,z.  " + A cos  * 


where 


A"  « -2 


r 

k 2E*  + 


c k ‘L  + k 1 
e 0 e e 


(z  +Z.) 

e 0 0 Jj  (xr)  dx 


(3.48) 


(3.49) 


The  transverse  electric  field  Intensities  In  the  air  are  determined  by  su' 
stltutlng  equations  (3- 46)  and  (3.48)  Into  equations  (3.23)-(3.28) . More 


specif  leal ly. 


E ® Id*  [-2.  (r)  + 


Z0 (r) 3 


(3.50) 


CHAPTER  4 

Coupling  Model  (Approximate  Solution) 

Our  next  task  is  to  evaluate  the  current  Induced  on  the  conductor  by 
the  lightning  strokes.  Practically  all  burled  cables  presently  In  use  have 
some  sort  of  plastic  jacket.  Even  if  a bare  conductor  Is  used,  there  will 
be  a layer  of  oxidation  at  the  surface  of  the  metal  wire.  In  addition,  It  Is 
likely  that  a thin  layer  of  air  Is  present  between  the  metal  surface  (or 
metal  oxide  surface)  and  the  soil.  In  short,  we  should  not  expect  the  bare 
wire  to  be  in  electrical  contact  with  the  soil.  To  be  realistic,  a lossy 
dielectric  layer  is  Included  in  our  consideration. 

Most  cables  are  buried  30"-36"  below  the  surface  of  the  earth.  For  most 
frequencies  of  practical  interest,  the  current  which  arrives  at  the  far  end 
of  the  wire  due  to  a lightning  stroke  is  essentially  insensitive  to  the  var- 
iation of  the  depth  of  the  cable,  provided  that  the  wire  is  long  and  that  the 
depth  is  deep  enough.  Thus,  In  this  section,  the  model  depicted  in  Figure  2b, 
which  is  an  approximation  to  the  configuration  shown  In  Figure  2a,  will  be  used 
to  derive  a set  of  differential  equations  governing  the  current  induced  on  the 
buried  wire. 

4.1  Transmission  Line  Equations 

The  axial  direction  of  the  wire  is  taken  to  be  the  x-axis,  (Figure  3). 
f he  field  at  any  arbitrary  point  in  the  space  can  be  decomposed  into  two 
parts.  One  part  is  the  field  produced  by  the  external  sources,  in  the  absence 
of  the  wire  and  the  insulation  layer,  and  will  be  designated  by  a superscript  0; 
the  second  part  is  due  to  the  current  induced  on  the  conducting  wire,  and  will 
be  designated  by  superscript  1. 

To  find  the  field  produced  by  the  current  on  the  wire,  we  start  with  one 
of  the  Maxwell's  equations  (2.1 ')  in  Integral  form: 


0 


dfi,  = -Ju 


dt 


(4.1) 


arid  choose  the  contour  to  be  C ^ , as  shown  in  Figure  3a.  As  shown  in  Figure 
3a,  the  line  integral  along  the  path  can  be  broken  into  five  parts 
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INSULATION 

€i  »Mi 

CONDUCTOR 
«c  »Mc  »°“c 


C L 

t 1 (r,x)dr  + E ' (r ,x)dr  + 
re  r i 

c 


/•x+Ax 


Ex,  (b,x)dx 


, 1 (r,x+Ax)dr  + jl (r,x+Ax)dr 
■*  c 


re 


= -jw 


x+Ax 


n: 


Bqj  (r ,x)dr  + 


BQ  1 (r ,x)drj  dx 


J 


After  collecting  the  terms  and  noting  that 


c r°°  i r®0  3E  * (r,x) 

E_  (r,x+Ax)dr  - Ax  — ^ dr 

< r 


l 


j Ere  <r’x)dr  + 


] 


‘ c 
c 


re 


Sx 


E (r,x)dr  + 
r i 


1 ~ rc3Erilv 

E^.  ( r,  x+Ax)dr  - Ax  j £ 


ri<r,x) 


dr 


we  have 
Ax 


(C  9Er|  (r,x) 


f°°3E (r,x) 


3x 


dr  + 


re 


3x 


dr  + E . (b 
xl 


•x)l 


= -jtoAx 


B-.  (r ,x)dr  + 


01 


6e 


1 (r,x)dr^J 


If  a voltage  V and  a flux  $ are  Introduced  as 


V (x)  * J Er 1 1 (r’x)dr  + 


1 


E (r,x)dr 
re 


0(x)  * Bol  (r ,x)dr  + 


c 

r°°„  1 


61 


V (r-x)dr 


then  the  equation  (A. 6)  becomes 


dV  (x) 


- j C0<|)  (x)  - E . (b,x) 


dX  J"T'"'  “2  I 

This  is  one  of  the  desired  differential  equations. 
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(A. 2) 


(4.3) 


(4.4) 

(4.5) 


(4.6) 


To  find  the  second  differential  equation,  we  examine  another  one  of 
the  Maxwell's  equations  (2. 2')  in  Integral  form: 


& H'*d&  ® ( °,+Jwe  )tf^  *d s 

Jr  J<;  > ‘ ! 


(^*7) 


and  choose  the  contour  to  be  C^,  as  shown  In  Figure  3b.  Clearly,  the  line 

Integral  along  path  AB  and  along  path  DE  cancel  exactly.  When  the  wire  Is  made 

of  good  conductors,  l.e.,  a »we  , the  displacement  current  term  In  the  wire 

c c 

may  be  neglected,  then 


•BCD^,d^=  'x(x+Ax) 

LAff5,dt  - “'x(x) 


(^.8) 


(M) 


where  I Is  the  axial  current  carried  by  the  wire.  Substituting  these  two 

X 

equations  Into  (A.1C),  we  have 
dl  (x)  , 

- ~(ai+jwc,)Er(  (b ,x)  * 27Tb  (A.lo) 

To  make  (A. 9)  and  (A. 11)  useful,  It  Is  necessary  to  express  V,  $,  E \ 

and  E ^ In  terms  of  I explicitly.  For  this  purpose,  three  Impedances  are 

, , * ,(8) 

Introduced  : 


(1)  Surface  Impedance  (or  skin  effect  Impedance) 
E (b ,x)  E °(b,x)  + Ej(b,x) 


X X 


(11.11 ) 


(II)  Longitudinal  (Inductive)  Impedance 

zl  - +J“r^r 

x 

(111)  Transverse  Impedance 

l e . 

T 2Trb(0,  tj(D£  j)Er  ! (b,x) 


(A.  12) 


('♦.13) 


Upon  substituting  (4. 11 )- (4. 13 ) Into  (4.6)  and  (4.9),  these  two 


equations  assume  the  form  of  the  transmission  line  equations. 


c'V(x) 

dx 


*'x(x)  (ZL  + V 4E0  {X) 


(4.14) 


(x) 


dx 


4 V(x) 

n 


(4.15) 


Figure  4 gives  a stmpie  and  familiar  representation  for  the  transmission  line 
equations.  As  we  are  mainly  interested  in  the  current  1^,  it  would  be  desirable 
to  eliminate  the  voltage  V in  favor  of  the  current  1^.  However,  the  resulting 
equation  would  be  extremely  complicated  as  the  Impedances  Z^,  Zs>  and  Zj  are  also 
dependent  cn  x.  Fortunately,  the  impedances  Z,.,  , and  Zj  are  approximately 

constants  under  certain  circumstances,  to  be  discussed  in  Appendix  A.  Under  the 
approximation  indicated  there,  these  Impedances  are 


ZL  - ZU  + ZL2 


ZT  * ZT1  * ZT2 


(4.16) 

(4.17) 


The  expressions  for  7.L]  , Z^,  Zj^  , Zj2 , and  Zg  are  given  In  Equation  (A.32)- 
(A.36)  While  these  equations  are  approximate,  they  are  valid  for  most  situations 
of  practical  interest.  In  addition,  from  the  definition  of,  and  the  expressions 
for,  the  longitudinal  and  transverse  impedance  Z(  and  Zj,  it  is  clear  that  these 
expressions  can  be  extended  to  include  multi-layer  cobles.  When  the  impedances 
Zs,  ZL,  and  ZT  are  approximately  constant,  equations  (4.14)  and  (4.15)  reduce 
to  the  well-known  transmission  line  equations  with  a distributed  source  term 

E (x).  The  solution  to  such  a system  of  coupled  differential  equations  can 

0 

best  be  expressed  Jn  terms  of  voltage  and  current  Green's  functions  as 
discussed  in  Chapter  6. 
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4.2  Distributed  Fields  Due  to  Lightning  Strokes  to  Ground 

We  shall  use  the  results  just  obtained  In  Chapter  3 to  calculate  the 
fields  and  the  potentials  of  the  lightning  stroke  to  ground.  The  lightning 
strokes  are  represented  by  a current  I along  the  2-axis.  The  current 
path  Is  subdivided  Into  small  segments  dz  and  the  Hertz  vector  tt  due  to  each 
segment  Is  given  by  equations  (3.4)  and  (3.12).  Thus,  the  total  Hertz 
vector,  for  z < 0,  Is 


TT 

ze 


r 

dz  1 

J« 


r 2ko 2t 
1 • ke\  + kA 


eJg (tr)  dr 


(4.18) 


and  the  scalar  potential  4>  Is,  by  substituting  the  above  equation  Into  (2.16), 


4>e(r,z) 


• 

J 2irKe 


r L 


zt 


0 J 


e 

X 


•°  'k.  50  + W 


j— TO  (tr)  dr 


(4,19) 


Of  particular  Interest  to  us  Is  the  case  of  low  frequencies  where 
toe  « a,  and  the  displacement  current  term  Is  negligible.  Under  such 
conditions,  a rather  simple  approximation  for  (4.19)  can  be  obtained.  By 
letting  toe-K),  we  obtain, 

*e<r'z)  ‘ 2no~  I eZT  J0(Tr)  dT 
e } o 

(7) 

The  integral  can  be  evaluated  and  leads  to 


<!*  Jr.z)  - 


2770 


[r 


TP71 


(4.20) 


which  is  the  same  as  the  potential  at  (r,z)  Inside  or  on  the  surface  of  the 
earth  when  a dc  current  I Is  entering  at  the  origin. 

In  particular,  the  distributed  source  Eq(x)  Is  given  by 
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CHAPTER  5 

Coupling  Model  (Exact  Solution) 


In  order  to  establish  exact  formulas  for  the  propagation  characteristics 
on  earth  return  conductors,  rigorous  theoretical  solutions  for  the  propagation 
of  current  along  an  extended  conductor  In  Infinite  contact  with  the  earth  are 
obtained  for  the  case  In  which  the  current  enters  the  earth  at  an  electrode 
which  represents  the  terminal  effects  of  a lightning  channel.  The  general 
case  Is  again  considered  here,  In  which  the  conductor  Is  not  necessarily  In 
direct  contact  wlththe  earth,  but  has  a contact  Impedance  wl th  the  earth, 
as  In  cables  with  coverings  provided  for  corrosion  or  mechanical  protection 
or  electrical  shields, 

5.1  Transmission  Line  Equations'. 

Consider  an  extended  straight  conductor  of  radius  p half  burled  In  the  xy 

o 

plane  of  a rectangular  coordinate  system  (x,y,z),  as  shown  In  Figure  5,  and  let 

the  x axis  extend  along  the  conductor.  Let  $ (x)  and  <f>  (x,y)  be  the  scalar 

c e 

potential  In  the  conductor  and  In  the  earth  (at  the  separation  distance  y) , 
respectively.  Also,  let  Ac(x)  and  Ag(x,y)  be  the  x components  of  the  vector 
potentials  In  the  conductor  and  In  the  earth  (at  the  separation  distance  y) . 
respectively.  It  Is  assumed  that  there  are  no  radial  variations  In  the  field 
variables  inside  the  conductor. 

The  x component  of  the  electric  field  Intensity  E^(x)  along  the  surface 


of  the  conductor  Is 

Ec(x)  ' • ♦cW 


(5.1) 


and  the  x component  of  the  electric  field  Intensity  Ee(x,Pg)  In  the  earth 
adjacent  to  the  conductor  surface  at  y^Pg  Is 


Ee<*’p0) 


dx 


<}>  (x,pn)  -jwA  (x,  Pn) 


(5.2) 


Since  there  Is  no  change  In  the  magnetic  flux  between  the  conductor  and  the 
earth  adjacent  to  It 


Ac(x)  » Ag(x,P0) 
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(5.3) 


i:rr\rr>t 


The  electric  field  intensity  Ec(x)  along  the  surface  of  the  conductor  Is  also 
given  by 

Ec(x)  - Zj  I c(x)  (5.4) 

where  Zj  Is  the  internal  Impedance  of  the  conductor  per  unit  length.  As  shown 
in  Appendix  B,  k } JQ(kP0) 

Zi  Jwe  2iiPQ  J,  (kpp) 

From  the  above  equations 

Ee*X,p0*~2i  'c^  5 dx  ■ (5.5) 

The  resultant  electric  field  Intensities  are  regarded  as  the  sum  of  an 

Impressed  primary  field  (denoted  by  a single  prime  superscript),  due  to  the 

lightning  channel,  and  an  induced  secondary  field  (denoted  by  a double  prime 

superscript),  due  to  current  In  the  conductor,  i.e.,  let 

E (x,y)  » E ‘ (x, y)  + E "(x,y)  (5.6) 

6 C C 

and 

<t>  (x,y)  * <t> 1 (x,y)  + 4>."(x,y)  (5.7) 

c c e 

<J>  (x)  “ <J>  1 (x)  + <P  "(x)  (5.8) 

c c c 


The  Induced  secondary  potential  between  the  conductor  and  an  adjacent  point  in  the 
earth  has  the  following  relation  to  the  leakage  current  t^(x)  where 

4>C'U)  -'4>e"(x,P0)  - l£(x)/  Y , (5.9) 

and 


lo  (x) 


dx 


'c(x) 


(5.10) 


and  Y I s the  admi  t tance  of  the  conductor  Insulation  per  unit  length.  Also,  for  a 
wire  of  infinite  length,  the  Induced  secondary  electric  field  Intensity  due 


to  a current  I (x)  distributed  along  the  entire  length  of  the  conductor 
c 


Is 


+00 


ee"<*.P0> 


«*> 


'c(x')+57 


dx,zt(rp„x)  a?1  'c(x,)  (5-n) 


l‘\T"  ■■■-— ** 


where 


_ / 2 , 2 

V p°  x 


(5.12) 


and 


-X  5 x-x* 


(5.13) 


and  2 0(r  ) Is  the  longitudinal  electric  field  Intensity  at  the  surface  of  the 

1 P0X 


conductor  at  x,  for  a unit  axial  current  at  x' , while  Z (r  ) Is  the  trans- 

t P0X 

verse  electric  field  intensity  at  the  surface  of  the  conductor  at  x,  for  a 
unit  radial  current  leaving  the  conductor  at  x*. 

Therefore,  by  superposition,  the  following  Integro-dl fferentlal  equation 


is  obtained  for  the  current  I : 

c 


d_ 

dx 


7.  sr  'c(x)  + 


-j-y  ^ 

dx'zt(V)  37-  lc(x',J  - zl'c 


(x) 


+00 


d*'V'Vx>  'c'*'1  ■ 'E0(x) 


(5.U) 


where  the  impressed  electric  field  Intensity  Eq(x)  along  the  conductor  Is 
related  to  the  impressed  electric  field  Intensity  E'(x,p.)  In  the  earth  next 

© v 

to  the  conductor  by 


E^x)  - - ^<Kx)  - e;(x,p0)  - ~ r^(x)  -^(x,P0)] 


(5.15) 


Equation  (5>14) Is  an  exact  generalization  of  the  approximate  transmission 

line  equations  (4,14)  and  (4.15)  derived  earlier.  Once  the  current  I (x) 

c 

is  determined  from  the  above  transmission  line  equation,  the  conductor 
potential  <j>  (x)  and  the  potential  In  the  earth  <f>  (x,y)  are  also  determined  as 

C G 

Integrals  over  the  current.  The  difference  between  the  Impressed  potential 
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4>0(x,y)  and  the  resultant  potential  <J>e(x,y)  in  the  earth,  which  Is  the 


negative  induced  secondary  earth  potential  <J>M(x,y),  is  given  by 

e 


(10) 


<f>(x,y)  - <P  (x,y)  = + 
0 c 


t +» 

dx 1 


•7  (r  )£_  | (x> ) 

V yx'dx'  V ' 


(5.16) 


Also,  the  induced  secondary  potential  4>" (x)  is  greater  than  the  induced 

c 

secondary  potential  4>y(x,y)  of  an  adjacent  point  in  the  earth  by  an  amount 


I (x)/Y  such  that 

J6  I 


(10) 


Vx’po}  “ *c(x)  e y: 


d_ 

dx 


t +oo 


I (x)  + 
c 


dxZt  (r  ) b 'c(x,) 


P0X 


(5.17) 


5.2  Integral  Transform  Representations 

In  order  to  solve  the  characteristic  equation  (5.14)  for  the  current  I (x) , 

c 

a Fourier  Integral  transformation  between  the  space  variable  x and  the  wave  number 
variable  ? is  Introduced,  i.e.,  let 


dxf{  x)e“J;x 


(5.18) 


(5.19) 


where  f and  ? are  any  of  the  field  variables  or  sources  defined  above.  Therefore, 

In  the  "wave  number"  domain,  the  characteristic  equation  for  the  transformed 

current  t (t;)  is 
c 
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f-  * cS 

‘V1  + zi  * «c<0  - •„(« 

(5.20) 

where 

- 

*t(V  * J 

W ml 

dx  h\l**SiX 

m w 

(5.21) 

* 

and 

fdxZ,(r  )e+JCx 

L p*x 

(5.22) 

'°W  ■ i ; 

dx  E0(x)e^x 

-00 

(5.23) 

Therefore 

1c(0  - e0(C)  / A(c) 

(5.2*0 

where 

r2 

A(0  •*  •^y  * 2 

(5.25) 

and 

h (, 

rP.C} 

(5.26) 

i + **(r0  r)  (5.27) 

Tl  c p0^ 

Therefore,  to  summarize  the  preceding  derivations,  the  current  I (x ) , the 

c 

conductor  potential  4>  (x),  and  the  potential  In  the  earth  4>(x,y)  are  given  by 

C fs 


, (x)  „ JL 
c fa 


•o.(c)  +JCx 


1 * 


W 


, r<°  en  (c) 

Mx)  ■ ~ ~i  j dc  5~EIc7 


fa 


*00 

00 


\T  + z*  (r«  r } 

Y,  t PQC 


+JCx 


I f " ®rt (?) 

Vm)  • Vx.y)  - ^ J zt(rP0c)e 


+JCX 


(5.28) 

(5.29) 

(5.30) 


5,3  Propagation  Constant 

To  find  the  current  I (x) , the  conductor  potential  4>  (x) , and  the  earth 

c c 

potential  <t>  (x,y) , the  following  transformations  of  the  transfer  and  longitudinal 

€5 

impedances  are  required. 

' i 


(5-51) 


zjr  *> 
t P05 


»+  CO 


dXZ  (r  )e+J?X 

1 PqX 


z Jr  ) 

^ Pq^ 


. +00 

r*  (v,e 


where 


/ 2 , 2 
/Da  + C 


rP05B/Po 

^(rp  j.)  and  zt(rP  are  is, 

respectively,  of  t$o  colli 

When  the  earth  Is  uni  ted, 

the  functions  z.(r  ) anc 
* PqX 

horizontal  dipole,  are 


Zt(rP0X)  ~ Jwe 

I e“YrPnX 

x f'Q*  — ■ p V 


whe  re 


Yrn.v  . .. 


Yrp 


y = jk  = [-u)2ue]1//2  » [jwp  (o+jwe)} 

Therefore,  the  functlonszf rQ  r)>'zl  (rn  t)  are 

r 0 0 


Z<V 


2n<rp0‘:>  " 2 


7T  Jans  VUIpo* 


J?|K,(k|P0)  - /TV 


Since  the  Important  part  of  the  Integration  range  z_(r  ) and  z (r  r) , as 

l PqC  t PqC 

given  above,  vary  nearly  logarithmically  with  £,  large  variations  In  £ produce 
small  changes  in  these  functions.  For  this  reason  it  is  permissible  to 
approximate  these  functions  at  C » where  Is  a constant  so  chosen  that, 
in  the  important  part  of  the  integration  range, 
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(5.38) 


**  V : z‘(rvo) 


z„(r  “ z/>(r„  r ) 

1 V £ P0C0 


(5.39) 


With  sufficient  accuracy  for  most  practical  applications  Sq  can  be  taken  equal 
to  f,  the  propagation  constant  of  the  cable,  In  which  case  the  following 
transcendental  equatton  for  T Is  obtained: 


I*  1) 

T 


» Z 


z~r 


c- r 


(5. AO) 


or 


r2  [I  + 1 

1 L V TT 


and 


r - 


1 

" jwe 

Ko< 

l<^2 

+ r2: 

- 1 

l 

TT 

JUK? 

~ 1 
2TT 

Jwy 

Jwy 

‘ “2tT 

In 

, Jwil 

z,  +r  — t|r|K  (|r|P„) 

i Y p 

0 


/y2  + r2  Kj  (/y2  + r2 


In 


1.12... 


In 


1.85... 


1.85... 

p,/y  2 V ? 


1 + L_  ,n 

VL  TTjuie  YpQ 


(5.A1) 


(5.A2) 


(5. A3) 


(5.AA) 


An  approximate  solution  to  the  above  transcendental  equation  Is,  by  Inspection, 
2 


p 2 ^ 


(5.A5) 


Inserting  this  value  for  Tlnto  the  logarithmic  terms  above,  the  expression 
reduces  to 


2 


TpT 


■r 


In 


(5.46) 


Since  the  ratio  of  the  above  logarithmic  terms  Is  practically  unity,  for  all 
practical  purposes 

r'?T  (5.47) 

5.4  Ground  Strokes 

The  propagation  characteristics  of  the  current  1 (x)  along  the  surface  of 

c 

a buried  cable  due  to  lightning  strokes  to  ground  In  the  vicinity  of  the  cable 
are  now  determined. 

5.4.1  Direct  Strike  (Arcing) 

A lightning  stroke  to  ground  may  arc  directly  to  a burled  cable  in  the 

vicinity  of  the  base  of  the  lightning  channel,  in  which  case,  virtually  all 

of  the  current  will  enter  the  sheath  near  the  stroke  point. 

When  a lightning  current  I enters  the  sheath  at  the  point  x *»  0 and  the 

s 

sheath  is  assumed  to  extend  indefinitely  in  opposite  directions  from  this  point, 

at 

r}x  | 


the  sheath  current  I (x)  at  the  distance  x is  given  approximately  by 

c 


I 


«c(x) 


(5.48) 


where  Tis  the  propagation  constant  of  the  earth/sheath  circuit,  as  described 
above. 

The  induced  electric  field  EQ(x)  along  the  inner  surface  of  the  sheath  is 


E0U)  -f-  z,  • 


- r|x| 


(5.49) 

where  Zj  Is  the  internal  surface  Impedance  of  the  conductors  per  unit  length, 
as  developed  in  Appendix  B. 

5.4,2  Indirect  Strike  (Conductive  Energization) 

Alternatively,  a lightning  stroke  to  ground  may  be  too  distant  from 
a buried  cable  in  the  vicinity  of  the  base  of  the  lightning  channel  to 
arc.  directly  to  the  cable;  however,  in  this  case  the  cable  can  still  be 
conductivi ly  energized  by  the  current  entering  the  ground  at  the  base 
of  the  lightning  discharge,  which  is  easily  represented  mathematical ly  by  an 
electrode  near  the  wire. 
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- . ’ ™ ».  — *.  „rth  , a.actroda  „ th.  separation  d.stance 
ps  conductor,  th.  pressed  secondary  earth  potent..,  a.onp  the 

conductor  ts  given  by 


<*>0  <X>P0)  * 2 (r  ) | 

, 1 PsX  s 

where 

V ^7*7 


(5.50) 


‘X  “ x-x 


Therefore,  the  Impressed  electric  field  Intensity  along  the  conductor  Is  then 
given  by 


E0  ^ 


Therefore, 


d 

Sr*o(W  * / IS-f- -^r  JCtt(r  )e‘J'x’le45I<  (5.51) 

— /35f  L 1 0*c  I 


P0(c) 


-JCzJr  (e'1^' 
C 


(5.52) 


The  conductor  current  f j 

« ’ the  C°"ductor  Potential  * , and  the  earth 

potential  *#(xfy)  are  then  gfven  by 


'c(x) 


(r  ) 

t p,  r 


K —SiL+Jz* 


^2"  j A(?) 


(5.53) 


4»c(x) 


I -is 

^ /2n  Hz) 


(5.5^) 


vfRn»'7i 


^U,y)  ‘ i I ^ — 'jL^\*U'lrn}  .*}<* 


(5.55) 


where 


P 5 + y 


(5.58) 


Therefore,  after  substitute  the  values  for  the  transformed  longitude, 

Impedance  ^(r^)  and  the  transformed  transverse  Impedance  z (r  ),  ]nt0 

equations  (5.53, -(5.55)  the  conductor  current  I (x>  rh„  a'  ^ ' 

urrent  icU/,  the  conductor  potential  <J>  (x) 

and  the  earth  potential  ♦.(x.y)  are  approximated  by 


,c(x)  = 


^ jug  ^ 


(5.56) 


*c(x)  r 1ST 


(5.57) 


4>e(*»y)  ■ 4>n(x,y)  — ~ r-j-4 


» rJ^-  T9—  -♦<% rp£)] 


(5.58) 


Vo™ 

"Jue  + Y(  Kq< 


where 


dl  — ■e't 
/t^  + 


The  function  4>  is  related  to  the  Bessel  and  Neumann  functions  as  follows 


•{10) 


<j>(u,v)  - Jo(v)1n~>-  ~Y„(v)  + I At 


(5.  59) 


where 


and 


and 


W a y'u2-f  v2 

A0  - 0 
Aj  *B  w 

l“l 


wu  -(l-l)v2At_2 


(\  - 2, 3, A,. ..) 


■\  x 

where  J0  and  Y0  denote,  respectively,  the  zero  order  Bessel  and  Neumann  functions 
of  the  complex  argument  v. 

A1  so 


Ko(v)  ■ ) o (v)  1 n (y)'  + £ 


09  /V \iX 


(fY 


t-o  Mr 


♦ (X) 


(5.60) 


and 


l0(v)  * Jo(lv) 

where  to  and  Ko  are  the  modified  Bessel  and  Neumann  functions,  and  iJj  ts  the 

di gamma  function. 
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CHAPTER  6 

Distributed  Source 


The  voltage  and  current  waves  on  the  outermost  conductor  of  a coaxial 
transmission  line  excited  by  a distributed  voltage  source  due  to  a lightning 
discharge  as  determined  either  approximately  (Chapter  4)  or  exactly  (Chapter  5) 
are  now  determined. 

6 . i Telegrapher's  Equation 

The  transmission  line  extends  for  a length  i along  the  x-axis  of  a rect- 
angular coordinate  system,  as  shown  in  Figure  6.  The  transmission  line  has  a 
terminating  Impedance  Z at  x“0  and  a terminating  Impedance  Z+  at  x«°£.  The 
transmission  line  is  excited  by  a distributed  voltage  source  V(x'). 

The  voltage  V(x)  and  the  current  I (x)  on  the  transmission  line  satisfy  the 
coupled  wave  equations,  c.f,  equations  (4.14)  and  (4 . 1 5 ) » 

1-  V = - Zl  + V (6.1) 

I » - YV  + I (6.2) 

dx 

where 


and 


Y = l/ZT 

V - E (x) 
I * 0° 


These  coupled  wave  equations  are  easily  solved  with  the  use  of  Green's  Functions 
for  the  voltage  and  current. 


J? Green's  Funct  1 o 


The  voltage  V(x)  and  the  current  I (x)  on  the  transmission  line  due  to  a" 


distributed  voltage  source  V(x')  are  determined  by  the  superposition  Integrals 


V(x) 

Kx) 


dx'G  (x,x')i/(x‘) 

0 V 

1 

dx'G.  (x,x')t/(x‘) 
0 • 


where  the  voltage  Green's  Function  G^(x,x')  and  the  current  Green's  Function 


G|(x,x')  satisfy  the  coupled  wave  equations 
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— 


(6.5) 


3T«,  - - 26,  ♦ «(»-»') 


« *1  ■ -VGv 


(6.6) 


where  6(x“x')  Is  the  Dirac  delta  function  at  x -x'. 

The  coupled  wave  equations  (6.5)  and  (6.6)  are  solved  simultaneously  to  yield 


Gv(x,x')  - \ 


V_(e+Yx  + rVYX) 


(x  < x') 


(6.7) 


V+(e"Yx  + rVYV2y*)  (x  > x') 


G.(x.x')  ml 


- f (e+YX  . fe“Yx) 

c 


(x  < x') 


(6.8  ) 


(.-YX  - r+e+Yx.-zyt) 

C 


(x  > x') 


where  the  propagation  constant  y and  the  characteristic  Impedance  2 are 

c 


defined  by 
Y “ JTl 
1 ■ 


and  the  reflection  coefficients  T~are  defined  by 
1 - Z 


r * 


z + z 


+ z+  " z 
n+  _ C 


Z+  + Z 


The  constants  V are  to  be  determined  from  the  boundary  conditions  which  relate  the 
continuity  of  the  current  at  x **  x'  and  the  discontinuity  of  the  voltage  at  x ■>  x ' , 


1 .e. , 
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(6.9) 

(6.10) 


Gj  (x«  + e,  x')  - G( (x-e,  x')  • 0 
Gv(x'  + e,  *')  - Gv(x«-e>  x')  - l 


The  boundary  conditions  yield  the  following  values  of  V : 

- rVY*V2Y* 


»+YX'  - rVY*' 


where 


A - - 2(i-r“rV2YJl) 


Therefore, 


Gv(x,x«) 


^ (e+YX  + r e YX) 

2(i-r  r e 2yZ) 


e+7x'  _ r“e"Y* 1 YX  + +yx  -2y£.. 
i-’+?vi> — 'e  + r e e ; 

2 (i - r r+e+2n 


(x  < X') 


(x  > x') 


e-YX'  , rfenx'e-2yA 

2i  (i-rrVSYi) 


rVYX 


) (x  < X*) 


G,(x,x‘) 


enX'  - rVY"'  (e-YX  _ r+e+Yxe-2yA)  (x  > x,j 


2z  (i~r“rV2Y£) 


(6.11) 

(6.12) 

(6.13) 


(6.U) 


CHAPTER  7 

Transfer  Functions 


Once  the  Induced  current  and  voltage  surges  on  the  outer  conductor  (p«p>) 
of  a coaxial  conductor  are  known  due  to  a nearby  lightning  discharge,  the  Induced 
current  and  voltage  surges  on  the  Inner  conductor  (p  - p<)  can  be  determined  via 
the  use  of  the  Impedance  transfer  functions  for  a coaxial  cable,  as  developed  In 
Appendix  B.  As  shown  In  the  Appendix  B, 


Ex  P“P>  Zee*ext  + Zel*lnt 

Ex^p«p<  Zle*ext  + Zll*lnt 


(7*0 

(7.2) 


where 


1 Aee 


zi 

z. 


and 


:rjwF 

|c 

A 

a 

N 

i 

1 

1 

-Jw£ 

2tt  p<p> 

1 

1 

1 

-jwE 

2ir  P<P> 

k 

1 

-Jw2 

2irp< 

1L 

A' 


and 


Aee  * MkP<)Jo(kP>)  + Ji(kp<jYo(kp>) 
A,,  ” Y1(kp>)J0(kp<)  + Jl(kp>)Yo(kp<) 

A - Jx  (kp>)Yj  (kp^,)  - Jx  (kP<)Yi  (kP>) 

A'  » Jx(kp<)Yl(kp>)  - Jx(kP>)Yx(kP<) 


If  the  coaxial  cable  Is  composed  of  several  concentric  shells  rather  than 
just  one,  as  depicted  above,  then  the  above  formulas  can  be  used  recursively  to 
obtain  the  current  and  voltages  on  the  Innermost  cylinder  due  to  the  current 
and  voltages  on  the  outermost  cylinder. 
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CHAPTER  8 
Conclusions 

The  above  theories  are  all  that  are  necessary  to  find  the  Induced  voltage 
and  current  surges  on  the  center  conductor  of  a coaxial  conductor  resulting 
from  the  penetration  of  the  outer  armour  and  the  Inner  shield  of  the  cable 
of  a transient  electro-magnetic  pulse  caused  by  a nearby  lightning  discharge 
to  earth. 

The  resulting  equations  are  being  programmed  for  numerical  solution  on  a 
digital  computer.  When  the  program  codes  are  completed  and  checked,  a parametric 
study  of  the  results  will  be  undertaken.  The  results  of  the  parametric  study 
will  be  presented  In  a forthcoming  report. 
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APPENDIX  A 


Surface,  Longitudinal  and  Transverse  Impedances 

In  this  Appendix,  the  Impedances  defined  In  (*4.11),  (4.12),  and  (4.13)  are 
evaluated  and  the  results  are  listed  In  (4. 19)- (4.23) . The  derivation  Is  put 
In  on  Appendix  so  that  It  would  not  Interrupt  the  continuity  and  the  development 
of  the  main  text.  In  addlton,  because  of  the  geometry  Involved  In  the  derivation, 
It  Is  convenient  to  take  a Fourier  transform  with  respect  to  x,  Instead  of  a Hankel 
transform  with  respect  to  r,  as  Is  the  case  In  the  main  text.  It  Is  felt  that 
the  least  confusion  Is  caused  by  removing  the  derivation  to  an  Appendix. 


A.j Hertz  Vector 

The  geometry  to  be  considered  Is  a conducting  wire  of  radius  b,  surrounded 
by  n coaxial  Insulating  layer  of  thickness  c-b  and  embedded  In  an  extended  region 
as  shown  In  Figure  A.I.  The  axis  of  the  wire  Is  taken  to  be  the  x-axls.  As 
demonstrated  later,  all  boundary  condltons  are  satisfied  by  using  the  x-component 
of  Hertz  vector  ir^,  alone.  In  terms  of  the  x component  of  the  Hertz  vector  it  , 
the  field  Intensities  t and  are 


t 


/s 

+ X 


a27r  . 

f X . 

( l — 


. 2 . 
k.  it  . ) 

i Xi 


(A.I) 


•> 

H . 

i 


where  tr 

x 


9it  . a 

+ — 41 
must  be  a solution  of 


(A. 2) 


V2iv  . + k.  2tt  . » 0 (A. 3) 

XI  i X l 

The  continuity  of  the  tangential  components  of  the  magnetic  field  Intensity,  H, 
at  r“b  and  r®c  requires  that  
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earth 


INSULATION 

LAYER 


Fig.  A.l.  Geometry  of  an  Insulated  Wire 


(°c  + lu£c'  ~ 


<°1  + lue,)  -jp- 


(ol  * toel>  “ST 


(o  + Iwe  ) ' 

e e 9r 


at  r“b 


at  r-c 


(A.*) 
(A. 5) 


In  the  above  equations,  the  subscript  I stands  for  the  quantities  or  terms  related 

to  the  Insulating  layer.  As  (A. If)  and  (A.5)  are  valid  for  all  values  of  x,  It  Is 

35tt  , 

obvious  that  (o,  + Iwe.)  — Qjrg^r  also  continuous  at  r«b  and  r»c,  which  Implies 

that  the  continuity  of  the  normal  components  of  the  conduction  and  the  displacement 

currents,  (a,  + Iwe.)  E . , Is  automatically  satisfied. 

■ t r t 

The  continuity  of  the  tangential  electric  field  Intensity  E at  r*=b  and  r“c 
requires  that 


+ k ztr  « 
c xc 


+ k.2TV  . 
1 xt 


at  r**b 


(A. 6) 


+ k227T 

.2  I Vl 


+ k 2n 
e xe 


at  r*c 


(A.  7) 


Since  the  Interfaces  r-b,  -wCx"^00,  and  r«c,  -®<x<°°,  are  the  entire 
cylindrical  surfaces,  It  Is  convenient  to  use  Fourier  transform  with  respect 
to  x, 


ff  (r,£,w)  - 

X nr~ 


TTx(r,x,w) 


1 

00 

-00 

1 

r 00 

ir  ( 

J * 

*.co 

,u>)e^*dx 


,w)e"J^xd£ 


(A. 8) 
(A. 9) 


Applying  the  Fourier  transform  to  (A. 3)  and  noting  that  it  Is  Independent  of  <t>, 
the  differential  equation  for  Is  greatly  simplified  to 

I 3(r  *i  ) +(k.2-52)^  . - 0 

r 3r  3r  i xi 

It  Is  Immediately  obvious  that  the  solutions  for  are  the  appropriate 
cylindrical  functions.  Thus,  for  the  conducting  wire,  (r<b) , the  Insulating 
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region  b<r<c,  and  the  earth  r>c,  the  appropriate  expressions  are,  respectively, 


t.  co 

* (r»x*w)  - *7=  F (^)J  (e  r)e"J?Xd5  , 


(r  < b) 


(A. JO) 


(r.x.w)  * 

Xl 


r 00 

[Fi(C)J  (£,r)  + F (C)Y  (5,r)}e"JCxdC,  (b<r<c)  (A. 11) 

JLoo  ° 1 2 0 1 


where 


it  (r,x,w)  - - 

xe 

F, . - A*  - V2." 


.00  (2  ) 

F3(C)H  (£r)e"JCxdC 

J-00  0 e 


(r>c) 


(A. 12) 


To  ensure  the  convergence  of  these  Integrals,  the  branch 


Re(O>0,  Ik  chosen.  F , F , F and  F are  yet  to  be  determined.  In  terms  of  these 

012  3 

functions,  the  axial  current  I carried  by  the  conducting  wire  can  be  written  as 


I (x,w)  * 2tt 


J (r,x,w) rdr 


“ /Zr  b (a  + i we  ) 
c c 


f E F (5)J  (K  b)e”JCxdC 


c o 1 c 


(A. 13) 


Mini 


v '-'(.a),  «&(x)  , E (b,x),  as  Introduced  In  Sec.  4.1,  can  also  be  expressed  In 


terms  or  r , F , F , and  F . Substituting  these  expressions  Into  (4.1 1)-(4.13 ) 

o i 2 3 


F (i)C2J  (fb)e'jCX 


c o c 


d? 


2nb(o  + jwe  ) | F (C)£  J (£  b)e“JCxdC 

C C J O Cl  c 


(A. 14) 


Z,  (x)  a 


-jwpj  (a j + JuiEj) 


(Fl(rj)[Jo(?|c)-Jo(4,b)]  + F2(c)[YoU|C)-YoU|C)l}e"J^Xd4 


- 1 


2ftb  (a  +Jwe  ) 
c c 


F (t)U  (E  b)e~jcxdc 

O C 1 c 


(A. 15) 


juip  (a  +Jwe  ) 
e e e 

2iTb(a  +Jcoe  ) 
c c 


/F3(C)Hn(2) (E  c)e"JCxdC 
o e 

/F  (C)E  JjU  b)e"JCxd^ 


(A. 16) 
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(A.  17) 


- / {FiUHJo^cJ-^Ujb)]  + F2(c)[Yo(5|c)  - Yo(C,c)]  }Ce”JCxdC 
2irb(a +Jwe  } F (^)JiU  b)e"^CxdC 

C C CO  c 


ft F3(c)H  (2)(^  c)e"KxdC 
Z (x)  05  ■ ~iL_ .. ..  fi,  — 

T2  2 irb  ( a + j we  ) £ F ( ? ) J ! ( C b ) e " -1  ^xd C (A. 18) 

c c c o c 

Strictly  speaking,  the  Impedances  are  functions  of  x,  as  exhibited  In  t\.l4)~(A.l8) . 
However,  only  a moderate  variation  with  respect  to  x over  the  major  portions  of 
the  wire  Is  expected,  except  near  the  points  where  1 varies  raptdly.  Furthermore, 

X 

It  can  be  shown  that  when  the  fields  under  consideration  are  travelling  waves, 
these  Impedances  are  truly  Independent  of  x.  Thus,  a good  approximation  to  the 
Impedances  can  be  obtained  by  considering  the  natural  modes  guided  by  the  coaxial 
configuration. 

A. 2 Natural  Modes 

For  the  natural  modes,  the  external  sources  or  excitations  are  absent  and  the 
fields  are  completely  specified  by  (A. 10)- (A. 12) . To  determine  F , F , F , F , 

0 12  3 

and  the  propagation  constant,  equations  (A. 10) - (A. 1 2) . are  required  to  satisfy  the 
boundary  conditions  (A.4)-(A,7).  The  resulting  expressions  are: 

(o  +jwe  )£  J i (£  b)F  (c)  ~ (o.+jwe.U,  [Jj  (^. b) F x ($)  + Yj(£ ,b)F2(c)l  * 0 (A . 1 9) 

c C C C O I 11  I 1 

(Oj+JioejJSjUitejC^U)  + Yj^cJFjU)]  - (ae+Jwee)^eH,(2)(^c)F3(i;)  - 0 (A. 20) 

VJo(Ccb)Fo(c)  " ^,2[d0(^,b)Fj(i;)  + Yo(^b)F2(c)j  - 0 (A. 21) 

t|aiJ0(^«c)F»(c)  + Yo(5,c)F2(c))  - £e2H0(2)(£ec)Fs(0  - 0 (A. 22) 

(A. 1 9)- (A. 22)  form  a set  of  homogeneous  algebraic  equations  for  F , F , F , F . 

0 12  9 

for  a set  of  nontrivial  solutions  to  exist,  It  Is  necessary  that 
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(A. 23) 


where 


A ■ A 
1 2 


A 
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A 
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a +jwe 
c c 

a(+Jwe( 

^c 

a +jaie 
o c 

*1 

^c 

o +ju>e 
e e 

Si 

Oj+jwej 

«e 

a +Juie 
e e 

*1 

a(+]we( 

r 

e 

J.(eob)Vo(eib)  - Jo(5cb)V?,b> 

H.(2)<5ec)Jo(5lc)  - Hc(2)(5ec)J.(5lc) 
Hl  (2>  (^ec)Y0(^lc)  " Ho(2)(5ec)Yiafc) 


(A. 24) 


(A. 25) 


Equation  (A. 23)  Is  the  dispersion  relation  and  Is  used  to  solve  for  the  propagation 

constant.  Also  obtained  are  the  dependence  of  F , F and  F on  F : 

12  3 0 


r J (0  0 **  b) 

r Tf)'  ~ oyPJSt,'  5,  j^b)  - W.U.b') 

(A. 26) 

Mr,!  -io  +i.a.  ) (,  Wi(5  b) 

c c c c 

Frn  -3  rr^-Joit,  5,  J i (£,b)  " Hi  (F,,b) 

(A. 27) 

Fjr.)  ..tjwrc  F(Ec  J,^b)[Jo(5lc)  - rYoU,c)] 

r-JO  " a(+ jut,  Ce2  Ho{2,(F,ec)[J1(^|b)  - IYjUjb)] 

(A. 28) 

whe  r • 

VTi>  Jx(F.  b)J  (5,b)  - J (E  b)JxU,b) 

a +j  uf.  t,  col  o c 1 

1 o +j,jc  r 

^•Trr-  r ' J’<V’>',o(5,b>  • ■'0V>ME,b> 

(A. 29) 
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Let  the  propagation  constant  determined  from  (A. 23)  be  denoted  by  then 
the  current  I on  the  wire  must  be  determined  by 


r(x,w) 


-Kqx 


(A. 30) 


where  I Is  a constant, 
o 

From  (A. 13)  and  (A. 30),  F can  now  be  determined, 

o 


Fo(0 


I «(?-«„) 
o o 


v^fT  b(o  +J«e  ) £ J i (£  b) 

c c c c 


(A. 31) 


The  evaluation  of  the  Integrals  In  (A. 1A) - (A. 18)  are  greatly  simplified  by  the 
presence  of  the  del  tat  function  In  (A. 31).  Thus, 


£ JnKb) 

c o c 


2-rrb(a  +ju>e  ) Jx(£  b) 
c c c 


(A. 32) 


-jtop, (Oj+Jwe,)  Fi(^)[Jo(£1c)-Jo(£|b)]  + Fa(c)CYo(5|c)-Yo(C|b)] 
2trb(a  +jwe  ) 


jcope(oe+Jwee)  F,(0 
2Trb(a  +J(oe  ) F^tO  V j j'fe  V) 


Vo(?)J^cb) 
H (2)(£  c) 


(A. 33) 


C-C. 


(A. 34) 


-1 

2Trb(o  +joje  ) 
c c 


Fi(c)[Jo(£Jc)>J0(£|b)]  + F2(c)[Yo(£|c)-Yo(£ib)] 


, MO  Ho(2)(£ec) 

2irb(a  +jwe  ) 7~(rJ  £ Ji(£  b) 

c c o c c 


£cFo(0M5cb) 


(A. 35) 


c-c. 


(A. 36) 


These  expressions  are  valid  In  general.  They  can  be  simplified  for  highly  conducting 
wires  by  replacing  the  cylindrical  functions  with  their  appropriate  approximations. 
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APPENDIX  B 


Internal  and  External  Impedances 

The  Internal  and  external  Impedances  of  cylindrical  conductors  are  required 
when  the  frequency  Is  such  that  the  skin  effect  Is  considered. 

For  a long  cylindrical  conductor  oriented  along  the  x-axts,  with  Internal 
or  external  coaxial  return,  the  Hertz  potential  if  Is  Independent  of  <j>  and  x and 
has  a component  In  the  x-dlrectlon  only.  Therefore,  in  equation  (2.14), 


n-0  (B.l) 

5-0  (B.2) 

Hence 

T •=  k (B.3) 

Therefore,  the  Hertz  potential  is  explicitly  chosen  to  be,  ir  « ir^x,  where 

ir  *»  c J (kp)  + c Y (kp)  (B.4) 

X 1 O 2 0 

where  c and  c are  arbitrary  constants  to  be  determined  from  the  boundary 
1 2 

conditions.  The  field  Intensities  are  determined  by 

E - k2TT  (B.5) 

^ X 

Hx  - 0 (B.6) 


and 


f - V.  TT 
t 3x  t x 


(B.7) 


t ^ ^ IT 

H « -Jue  kxV  tt 

L L X 

Therefore,  the  only  non-zero  components  of  the  field  Intensities  are 


(B.8) 


E - k2[c  J (kp)  + c Y (kp)]  (B.9) 

X 1 0 2 0 

* -Jwe  k[c  J (kp)  + c Y (kp)]  (B.IO) 

V 11  2 1 

If  the  conductor  has  an  inside  radius  p<  and  an  outside  radius  p>#  the  magnetic 
field  Intensity  at  p ■ p<  must  be  zero,  provided  the  conductor  current  lftXt 

returns  outside  of  the  conductor,  and  at  p^p,^  must  equal  l.e.,  from 

Ampere's  Law  ...  . — - — — - 
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(B.n) 


2nP<H*| 


2irp>H 


*1 


P“P, 


P“P> 


ext 


(B. 12) 


Hence,  with  the  values  of  H from  the  above  equation,  the  boundary  eqndlttons 


for  c and  c are 
i 2 


<t> 


c J (kp  ) - c Y (kp  ) * 0 
1 1 < 2 1 < 


c J (kp  ) - c Y (kp  ) « “".-"V  v 
i i > 2 i > jwe  k 2irp> 

which  are  easily  solved  for  c and  c to  yield 

1 2 


c » 
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, I Y (kp.) 
1 ext  i < 

jwe  k 2irp>  A 


(B.13) 
(B. 14) 


(B.1 5) 


c = 


i I . J (kp.) 
J ext  i r< 

z -Ju)e“  k 2Trp>  A 

where  the  determinant  A Is  defined  by 


(B.16) 


A=  (kp>)Y^  (kp,^)  - J ^ (kp<)Y^  (kp>) 
Therefore,  the  Hertz  vector  becomes 


I 


IT 


ext 


Y1(kp<)J0(kP)  * J1(kp<)Yo(kp) 


x juc  k 2np>  A 

and  the  electric  field  Intensity  In  the  x direction  becomes 


E »■  - 
x 


'r,<kP<>Vko)  ' J.<k0<>Vk,>> 


j bit  2-rrp> 


(B. 1 7) 


(B.18) 


(B.19) 


The  ratio  of  the  electric  field  along  the  outer  surface  p ■ P>  to  the  total 
external  current  is  defined  as  the  external  surface  impedance  with  external 
return  t^\  and. the  corresponding  ratio  of  the  electric  field  along  the  Inner 
surface  p « p<  to  the  total  external  current  Is  defined  as  the  Internal  surface 
impedance  with  external  return  Z,  , I»e. 

I 6 


58 


z “ — r 

ee  ! 


Z.  « 
ie 


x I p-ps 


x ! p"p< 


k_ 1 jet 

jai£  2irp>  A 


__1_  1 

‘Jwe  2tt  p p>  A 


where 


Aee  “ Y1(kP<)jo(kp>)  ' WW 
The  following  Wronsklan  relation  has  been  used: 


(B.20) 


(8.21) 


(B.22) 


- VkP<)VkP<>  + JI<kP<)Yo(kP<)  ■ (B .23) 

If  the  conductor  hasan  Inside  radius  p<  and  outside  radius  P>t  the  magnetic 

field  Intensity  at  pap>  must  be  zero,  provided  the  conductor  current  ljnt 

returns  Inside  of  the  conductor,  and  at  p»p  must  equal  I /2irp.,  I .e. , from 

< Int  < 


Ampere's  Law 


2irp>H 


1p“P> 


(B.24) 


2Trp<% 


(B.Z5) 


Hence,  with  the  values  of  H from  the  above  equation,  the  boundary  conditions 


for  c and  c are 
1 2 


c J(kp')  - c Y (kp  ) - 0 
11*  2 i > , | 

c J (kp  ) ~ c Y (kp  ) * — j~7rr 
i j < 2 i < jwE  k Zlrp< 

which  are  easily  solved  for  c,  and  c to  yield 

1 2 

i Sot  Mkp>> 

Cl  “ -JuS  k 2^  A' 

i ixnt  J-(kp>) 

c2  jcoe  k 2 np<  A< 


(B. 26) 
(B . 27) 


(B . 28) 


(B.29) 


where  the  determinant  A'  Is  defined  by 


A'  - J (kpjY  (kpJ-J  (kp.)Y  (kpj 
i < 1 * i > i < 

Therefore,  the  Hertz  vector  becomes 


(B. 30) 
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x jwe  k 2KP 


W V^V^-Mkp^kp) 

A' 


(B.31) 


and  the  electric  field  Intensity  In  the  x direction  becomes 


I 


. . k * int 

x -jtue 


(B.32) 


2trP<  A, 

The  ratio  of  the  electric  field  along  the  outer  surface  p»p>  to  the  total  Internal 
current  Is  defined  as  the  external  surface  Impedance  with  Internal  return  Z^;  and> 
the  corresponding  ratio  of  the  electric  field  along  the  Inner  surface  p«p<  to  the 
total  Internal  current  Is  defined  as  the  Internal  surface  Impedance  with  Internal 


.return  Z{  f , I .e.  , 
E. 

Z. 


x P-P, 
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-Jw£  2irp< 


A* 


(B. 33) 


‘x|p-p> 


"el  I 


Int 


jwe  2tt  p<p>  l, 


0.34) 


where 


A',l  - (kp>)Jo(kp<)-Ji  (kp>)Yo(kp<) 

The  fol  lowing  Wronsklan  relation  has  been  used: 
' Yl(kP>)^0(^>)+Jl(kp>)Yo(kp>)  - 


(B.35) 


(B. 36) 


Therefore,  by  superposition, 


E - Z I +Z  t I,  . 

x p»p.  ee  ext  el  Int 


: I ■ Z,  I +Z,  , I . . 
x' P«P<  le  ext  1 1 Int 


(B. 37) 
(B. 38) 


Notice  that 


,1 


Therefore 
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For  solid  conductors  (p<  = 0) , Zge  is  usually  referred  to  as  the  internal 
impedance  Z.  of  the  conductor,  i.e. 

Z.  = !im  Z = -Ji~  -I- 

' p<-0  ee  jW£c  27rp>  Ji(kp>) 

When  the  Bessel  and  Neumann  functions  are  replaced  by  their  asymptotic 
expansions  for  large  values  of  their  arguments,  the  following  approximate 
formulas  are  obtained,  which  apply  for  the  type  of  cylindrical  conductors 
ordinarily  encountered. 

(B.40) 

(B.41) 
(B.42) 

(B.43) 
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where 


t - P>-  P< 
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APPENDIX  C 

Self  and  Mutual  Impedances 


The  self  and  mutual  Impedances  of  earth  return  conductors  are  required  between 
two  current  carrying  paths,  such  as  the  paths  C and  C'as  shown  In  Figure  Cl. 

The  fields  due  to  a Hertzian  dipole  located  on  the  curve  C'  and  oriented 
In  the  direction  of  the  z'  axis  have  the  forr. 


E 1 
z 


" !&'  t”Zp  (R)  + — , Z <R)] 

1 3z  2 t 


(C.l) 


E 


t 


IV—  V Z (R) 
Sis'  t t 


(C.2) 


Therefore,  the  fields  due  to  a continuous  distribution  of  Inflnl tlslmal  dipoles 
distributed  along  the  curve  C'  have  the  form 


, b'  a2 

E - 1 f ctt.  • [-Z- (R)  + Z (R)] 

z J.  3z*2  1 

a 

c'b' 

E ’ - 1 / dfc*  -V  V Z (R) 

t a'  3z  * * 

C' 

The  voltage  V Impressed  along  the  curve  C In  the  same  plane  as  C* 


(C.3) 

(C.A) 

Is  determined  by 


b b 

V « - / c£  *E  * - / dz(z*z'E  ' + z*t  1 ) 
a a z t 

C C 

Therefore, 


(C.5) 


and 

3 3 , 3 , 

a!  " aF'  008'**  + air  9,n  ^ 

Tho  voltage  V Is  now  evaluated  as 


b b' 

V “ I / di.  / dJt  • 2.  (R)  cos<J>  -IZ^R) 
= - 1 * c 


a 

C 


a 

C' 


(b-a)  (b 1 —a  * ) 

The  mutual  Impedance  between  the  two  current  paths  c and  c‘  Is  defined  as 


The  double  Integral  Is  the  mutual  Impedance  between  the  current  paths  and 
represents  the  longitudinal  Impedance  between  the  wires,  while  the  second  te 
Is  Independent  of  the  current  paths  and  represents  the  transverse  Impedance 
between  the  wire  terminals  through  the  surrounding  medium. 
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and 


n 

The  vol 


8 8 
• - gp-  cos<}»  + sin  <}) 

tage  V is  now  evaluated  as 


V - I 
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/ 

a 
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b 1 

dft  / 6V  Z.  (R)  cos4>  - |Z  (R) 

a ' * c 

C' 


(b-a) (b  * -a  * ) 


The  mutual  impedance  between  the  two  current  paths  c and  c'  is  defined  as 


(C  * 7) 
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b b 1 

J dft  / dS,'  Z£(R)cos4>  - Z (R) 
a a 1 

C C‘ 


(b-a) (b*-a») 


(C.8) 


The  double  integral  is  the  mutual  impedance  between  the  current  paths  and 
represents  the  longitudinal  Impedance  between  the  wires,  while  the  second  term 
is  independent  of  the  current  paths  and  represents  the  transverse  impedance 
between  the  wire  terminals  through  the  surrounding  medium. 
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APPENDIX  D 


FAA  Lightning  Study  Participants 


Florida  Institute  of  Technology  - Cable  Testing 

Dr.  Andrew  W.  Revay,  Jr. 

Richard  M.  Cosel 


Lightning  and  Transient  Research  Institute  - High  Voltage  Facility 
James  C.  Stahmann 

Purdue  University  - Protective  Devices 

Warren  Peele  (Project  Leader) 

Dr.  Chln-Lln  Chen 


Georgia  Institute  of  Technology  - Equipment  Analyses 

Keith  Huddleston 
Dr.  Ronal  Larson 
Dr.  John  Nordgard 


Air  Force  Institute  of  Technology  - Reliability  Aspects 

Lt.  Col,  Jerry  L.  Hanson 
Professor  T.  L.  Regul Insk! 


